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We investigate the cliirality dependence of physical properties of nanotubes which are wrapped 
by the planar hexagonal lattice including graphite and boron nitride sheet, and reveal its symmetry 
origin. The observables under consideration are of scalar, vector and tensor types. These exact 
chirality dependence obtained are useful to verify the experimental and numerical results and propose 
accurate empirical formulas. Some important features of physical quantities can also be extracted 
by only considering the symmetry restrictions without complicated calculations. 
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Since the discovery of carbon nanotube(CNT) Jj, there 
have been extensive investigations on the unusual phys- 
ical properties of this novel nano-material. The sim- 
plest CNT is the single wall carbon nanotube (C-SWNT) 
consisting of only one rolled up graphite sheet, which 
was first synthesized in 1993y, yj and now can be 
produced in large scale 0, 0]. As an analogue to the 
graphite, the III-V layered compound, namely, boron 
nitride(BN) sheet has also the hexagonal lattice struc- 
ture and can be wrapped into various nanotubes too. 
The single wall BN nanotube(BN-SWNT) was synthe- 
sized in 19960, which also attracts much attentions very 
recentlyHQ. Unhke the nonpolar C-SWNT which could 
be metallic or semiconducting with a moderate gapQ, 
this heteropolar nanotube is expected to be always a 
wide gap semiconductor [lol . Both of these two kinds 
of nanotubes have the similar descriptions of their chiral 
structures and constitute promising materials for wide 
applications 0. 

The properties of nanotubes are determined by their 
chiral structures. Therefore, to investigate the chirality 
dependence of various physical quantities of nanotubes is 
an interesting topic all the time. Particularly, to analyze 
the chiral composition of bulk samples will be helpful to 
the production of nanotubes with different species. Nu- 
merous experimental and theoretical investi gati ons are 
devoted to this subject HI El El E 113, El 111 El 
li^ l . Most recently, it has been reported that detailed chi- 
rality distributions in the bulk samples of the separated 
C-SWNT can be obtained by the fluorescence measure- 
ment and by the resonant Raman spectroscopy^! 
Their assignments of the chiral number (rii , 712) to the ob- 
served spectra are based upon the comparison between 
the tight-binding (TB) calculations and the experimental 
data. 

In the fluorescence experiment E3, the authors also 
gave the empirical formulas of the measured von-Hove 
singularities with respect to the structures of C-SWNTs, 
which are expressed in terms of the chiral angle 9 and 



chiral index v defined as, 9 = atanA/37i2/(2ni -I- n^), and 
V = mod(ni — n2, 3). Actually, for the purpose to give the 
dependence of physical observables on the chiral struc- 
ture of nanotubes, it is more suitable and efficient to use 
v and 9 instead of the chiral numbers (ni,n2). Since 
there are many types of nanotubes not being observed 
in the experiments or numerical calculations, an appro- 
priate empirical formula are quite useful and convenient 
to predict the properties of those unobserved ones. Gen- 
erally speaking, there are two possible methods to ob- 
tain these empirical formulas for different physical quan- 
tities. One is through fitting the numerical and experi- 
mental data, the other is from the analytical expansion 
around the two Dirac points of the hexagonal Brillouin 
zone based on various TB models and effective mass ap- 
proximations. However it has technical difficulties in the 
second way to get higher order terms which are some- 
times important. 

In this paper, we show the symmetry restrictions on 
the general chirality dependence of physical quantities 
of various types such as scalar, vector and tensor. This 
leads to compact forms of the chirality dependence for 
these observables on 9 and v. Since the results are model 
independent and exact, they not only can be used to pro- 
pose accurate empirical formulas from the numerical or 
experimental data, but also can indicate some important 
features of physical quantities without complicated calcu- 
lations. This idea was originated in our previous study of 
the natural optical activity of C-SWNT |22|. Our symme- 
try analysis is essentially based on the hexagonal struc- 
ture, and we will consider its application to the proper- 
ties of both the heteropolar BN-SWNT with large ion- 
icity and nonpolar C-SWNT. We examine the following 
physical quantities as concrete examples to illustrate our 
method, such as the excitation gap, electric polarization, 
dielectric tensor, and piezoelectricity. 

Considering a hexagonal lattice with base vectors oi 
and a2(see Fig. (^l), it can be rolled up into a nanotube 
along the chiral vector R = nidi + ri2a2|l2|, so that each 
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nanotube can be simply represented by a pair of chiral 
numbers (ni,n2) and the chiral angle 9 is just the angle 
between R and Oi. It then can be established a map- 
ping / from the space of chiral vectors on the planar 
sheet to that of the nanotube structure in a fixed way of 
wrapping. However this mapping is not one-to-one. For 
a given chiral vector Rq with chiral angle on the BN 
sheet with C31, symmetry, there are another 2 equivalent 
vectors, obtained by rotating Rq by 27r/3 successively, 
corresponding to the same nanotube. These vectors, as 
shown in Fig. have the explicit forms 

^0 = nidi + 722(32, 

i?2 = «2ai - ("1 + ?32)a2, • (1) 
Ri = —{ni + n2)di + nid2 

All of them constitute an invariant subspace of the three 
fold symmetry of BN sheet with the same chiral index v. 
Another three chiral vectors R^ and i?5 in Fig. 
have chiral angle 0-|-7r/3, O + n and O + hn/i, respectively, 
and can be written as 

^1 = («1 + "-2)ai - "102, 

Rz^ -nidi-n2d2, ■ (2) 

i?5 = -'^201 + ("1 + "-2)52 

They also form an invariant subspace, but with an oppo- 
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FIG. 1: Illustration of hexagonal lattices. The solid and open 
circles represent boron and nitride atoms for the BN-SWNT 
and both are the carbon atoms for the C-SWNT. 

site chiral index —v. The nanotube mapping from Ri is 
related to that from i?o by rotating the latter tube upside 
down. In addition, there is another special chiral vector 
R'o = ("-1 + ^2)0,1 — ri2a2, which is the reflection of i?o 
about the base vector di on the BN sheet and has the 
same chiral index v oi Rq. When mapping onto a nan- 
otube, it corresponds to the mirror image of that from 
Rq with respect to the section along the tube axis. For 
the graphite sheet, since the two atoms in one unit cell 



are the very same, it has a higher symmetry C^^. The six 
vectors Ri{i = 0, • • • , 5) all correspond to the same nan- 
otube which means that a C-SWNT may be represented 
by opposite chiral indices, v and —v. 

Briefly, through the mapping / the manipulations on 
R in the chiral vector set lead to the change of the struc- 
ture in the nanotube set, and we have the following three 
observations 

i. when 9 — > 6+2TT/3, the nanotube keeps unchanged, 
and v is also unchanged. 

ii. when 9 9 + 7r/3, the nanotube is rotated upside 
down, and v ^ —v. 

iii. when 9 —9, the nanotube is reflected with re- 
spect to the section along the tube axis, and v is 
unchanged. 

According to the first observation, we know that the 
physical quantities of nanotube can always be expanded 
in terms of the triangle function of 9 with respect to each 
class of v 

00 

g('') (9) = J2 a^"^ cos(3n6') + b^i'^ sm{3n9) . (3) 

n=0 

Q is some physical quantity. The coefficients al'^'''s and 
fei^^'s are functions of those chirality independent vari- 
ables, such as the tube radius r and some external pa- 
rameters. Our analysis will reveal the characteristic role 
of the chiral index v in classifying the chirality depen- 
dence. In the following we will consider some examples 
to show that the second and third observations together 
with features of the physical quantity under considera- 
tion will reduce the above chiral expressions Eq.Q. For 
the convenience of discussion, the Cartesian coordinates 
for nanotube is introduced as: the tube axis is set as z 
direction and the cross section is the xy plane with x axis 
passing through one atom on the tube surface. 

1. Excitation Gap A. We can treat the band gap for 
the C-SWNT and BN-SWNT as a scalar since its value 
does not change under rotation or mirror reflection of the 
nanotube. Hence, its ^-dependence must satisfy 

AM(6' + 7r/3) = A''-'\9) 

A^''\-9) ^ A^''\9) . (4) 

Notice that in the first expression in Eq.Q, becomes 
— v when 9-^9 + tt/3. Then, from Eq.Q the chirality 
dependence of A reads 

A'-^\9) = aa±aicos{W) +02 008(69) + ■■■ 
A(0)(6l) = ao + ai cos(66') as cos(120) + • • • . (5) 

Clearly, the coefficients of cos(36') ior v — ±1 should have 
the same value but opposite signs. 
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For the semiconducting C-SWNT with v = ±1, the 
longitudinal optical excitation can be measured in the 
fluorescence experiment [TJ and is quite useful in analyz- 
ing the chiral composition of bulk samples. In Ref.[13j, 
the authors fitted the experimental data of the von-Hove 
singularities for v — \ and —1 separately. According to 
their fitting function, the absolute values of a\ are quite 
different for v = \ and v — which does not agree with 
the symmetry analysis Eq.© and suggests that it may 
need to consider higher order terms like cos (60) for more 
accurate fitting functions. In fact we have tried to fit all 
their data for both v = ±1 by just one four-parameters 
formula with a cos(66) term 
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with parameters p.^s to be determined, which satisfies the 
symmetry restriction Eq.lQ. It turns out that the result 
has a smaller root-mean-square deviation than that in 
Ref. 0. 

Recently, an interesting temperature dependence of 
bandgap A(T) was reported for the semiconducting C- 
SWNT in Ref.j^, that is when 9 is small, the tempera- 
ture dependence of gap is monotonic for v = 1 and non- 
monotonic for I' = ~1. The different behaviors of A(T) 
for ly — ±1 exist even for very close 9. This could be un- 
derstood by assuming that both and ai are decreasing 
functions of temperature. Then from Eq.Q, it is clear 
that the bandgap for v — 1 is monotonically decreasing 
with temperature. For — ^1 and small 9, the sign of oi 
is negative, so that A can be nonmonotonic as temper- 
ature varying as a result of the interplay between ao{T) 
and ai(T). However when 9 is close to tt/6, the cos(30) 
is nearly vanishing and only oq takes the dominant role, 
so that the temperature dependence of bandgap should 
behave similarly for both z/ = 1 and —1. 

2. Electric Polarization (EP) P . The macroscopic elec- 
tric polarization along the nanotube axis is the conse- 
quence of the broken sublattice symmetry of BN-SWNT, 
which was studied as a geometric phase in Ref. '25^ . They 
found the sign and size of the longitudinal polarization 
are dramatically dependent on the chiral structure of 
nanotube. This in fact has its symmetry origin, and we 
will show below that some remarkable properties of EP 
can be extracted by the symmetry analysis. 

Due to the helical symmetry of the nanotubejiH. l27l|. 
this EP vector only exists in the tube axis direction(z- 
axis) and vanishes in the cross section(a;?/ plane), i.e., 
Pz ^ and Px,y = 0. According to observations ii and 
iii and the vector nature of EP, we have 

Pi^) (0 + 7r/3) = ~Pi-''\9) 
P^''\~9) = PM((?) , 

which leads to the chirality dependence of Pz 

P^^\9) = ai cos(36') -I- ag cos(96l) • • ■ 



(6) 



Pi'^\9) = ±ao -I- ai cos(36') ± 02 008(66*) -(- • • ■ 



(7) 



Eq.(|Zj) shows evidently that the armchair tubes(ni — 
ri2 = N) with 9 = 7r/6 and v = Q have no elec- 
tric polarization. Another type of achiral tube, zigzag 
tubes(ni — N,n2 — 0), show different pictures. If we 
assume reasonably the coefficients of higher order terms 
are small enough comparing to the zero order term ao, we 
then obtain P^°'> w for mod(A^, 3) = 0, and P^^'' « ±ao 
for mod(A^, 3) = ±1, i.e., when N is increasing, the EP is 
oscillating among 1, and —1, which is just the striking 
result of BN-SWNT found in Ref.ji^. 

For the C-SWNT, the six fold Cqv symmetry in the 
chiral vector set guarantees Pz°\9 + 7r/3) = Pz~'^\9), 
which combined with Eq. © leads to Pi'^'' — 0, i.e., no 
EP in C-SWNT. 

3. Dielectric tensor e. This second rank tensor is re- 
stricted to have the following form by the helical symme- 
try of nanotubefirf: 



(8) 



with Cxx = eyy ai 
elements of e, denoted by en with i — x,y,z, have the 
same chirality dependence. Similar to the analysis of EP 
vector, we obtain 
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+ n/3)^e[r\9) 



ei-^{9 + n/i) = -e^-/\9) 



(9) 



by noticing that the diagonal term is unchanged when the 
tube is reversed or reflected, and the off-diagonal term get 
its sign changed. Then 

Ml 



flo ± fli COS(30) + 02 cos( 

'{9) = flo -I- 02 008(66*) 04 008(126*) 
bi sm{39) ± &2 8in(66') + • • • 
eW(5() = 6i sin(36l) -I- &3 8in(96l) • • • 



.(±) 

^xy 
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o.-:,K9) 



(10) 



The coefficients o^s and h'^s in the above four expressions 
have no direct relationships. Obviously, the diagonal and 
off-diagonal terms have quite different chirality depen- 
dence and the off-diagonal terms vanish for the zigzag 
tube whose chiral angle is 0. 

The discussion above is for the heteropolar BN-SWNT. 
For the C-SWNT, the higher symmetry requires e^J{9 + 
7r/3) — exy'^\9), which together with Eq. © leads to 
txy = for any kind of C-SWNT. 

4. Piezoelectricity e. Piezoelectricity is the response of 
the EP of the material to the mechanical strain, which 
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is a third rank tensor defined by the derivative of EP 
vector with respect to the elastic strain tensor u as 
ejjfc — dPi/dujk- For the quasi-one-dimensional nan- 
otube, we are concerned with the response of EP along 
the z direction to uniaxial(s) and torsional (t) strains, 

= dPjdu,, et = dP,/dut . (11) 

Us is the stretch strain along the tube, and Ut is the 
torsional strain around the tube circumference. They can 
be related to the second rank tensor Uij in the Cartesian 
coordinates system through equations Us — u^z and Ut = 
{xUzy — yuzx) /^'^ with = + y^. Then, it is clear that 
how Us and Ut transform under the rotation of the tube 
upside down or the mirror reflection with respect to the 
section along the tube axis. Consequently, we have 

ei^)(e + 7r/3) = -ei-^)(0) 

ei''\~9)^-ei''\e), (12) 

and the chirality dependence of piezoelectricity then 
reads 

e'-s^He) = ±ao -I- ai cos(36') ± 02 cos(66') + 

ef") {9) = ai cos{W) + as cos{9e) + ■■■ 

e[^\e) = 61 sin(36') ± &2 sin(66l) -1- 63 sin(96l) -I- • • • 

e|°^ (6*) = bi sm{39) + 63 sin(96') + • • • . (13) 

Eq. H13|l implies the zigzag(^? — 0) tube can only have 
piezoelectric response to the longitudinal stretch and no 
response to the torsion strain around the circumference. 
On the contrary, the armchair {9 = 7r/6 and v = 0) tube 
has response to torsion strain but not to stretch strain. 
These conclusions are in consistent with the numerical 
results of the ab initio and TB calculations in Ref.j^l, 
where the chiral angle by definition has a difference 7r/2 
from ours. Apart from this difference in the 9 defini- 
tion, in Ref.|2J| only cos(30) and sin(36') terms appear in 
the chirality dependence of Cs and et, respectively, while 
there are extra terms according to our results Eq. (|13|1 . 
In fact, their 0-dependence of Cs and et do not agree with 
their numerical data very well. The reason for this differ- 
ence is that in Ref. '2^ the chirality dependence of e* and 
Cs is inherited from the BN planar sheet [2^, so that only 
the W terms are permitted by the C^y symmetry of BN 
sheet. However when the sheet is rolled up, this planar 
symmetry is broken, therefore those C3„-forbidden 
terms are not forbidden any more, although they should 
be small by the continuous argument from sheet to nan- 
otube, i.e., they vanish when the tube radius tends to 
infinity. 

For the C-SWNT, the symmetry leads to additional re- 
striction e2^''(6' + 7r/3) = e^'^\9), which combined with 



Eq. (|12|l requires Ca (9) = for both a = s and t, 
namely, for C-SWNT there should be no piezoelectricity 
due to its nonpolar feature. 

In the above discussion, the expansion coefficients 
ai^^ 's and 's can not be determined by the symmetry 
argument which in fact depend on the chirality indepen- 
dent parameters, such as the tube diameter, magnetic 
flux, temperature, and so on. By tuning the external 
parameters, one can adjust the magnitude of ai^^^'s and 
61"'' 's, which will be helpful to identify the chirality of 
the tubes. As examples, one could refer to Ref. [2^ 
to find the different dependence of these coefficients on 
magnetic fiux or temperature. 

As a conclusion, we give the explicit ^-dependence of 
physical quantities for different values of by a symmetry 
analysis. It shows clearly that the chiral index 1^ plays a 
characteristic role in describing the chirality dependence. 
This model independent method may be used to verify 
the numerical and experimental data and also can give 
rise to some important properties qualitatively without 
complicated calculations. In addition, this method is not 
restricted to the examples illustrated in this paper and 
could be extended to other situations. 

The author F. Ye would like to thank Dr. J. Wu at 
CASTU for helpful discussion. 
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